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I. INTRODUCTION 

Observation of superconductor's response to a weak 
external electromagnetic stimulus allows us to collect a 
number of subtle characteristics 2 . The way the magnetic 
field is expelled from a superconducting (SC) sample in 
the spectacular Meissner effect can be used to infer about 
many fundamental features of the system. Therefore the 
phenomena at the length scale of the magnetic field pene- 
tration depth A, i.e. in the region at the edge of the sam- 
ple where the induced supercurrents effectively screen the 
external magnetic field, are subject to intensive studies 
both on the theoretical and the experimental fronts of 
the research in cuprate superconductors^. In particu- 
lar, the magnetic field penetration depth can be used as 
a probe of the pairing symmetry since it can distinguish 
between a fully gapped and a nodal quasiparticle exci- 
tation spectrum^ 4 -. The former results in the thermally 
activated (exponential) temperature dependence of the 
magnetic field penetration depth, whereas the latter one 
implies a power law behavior. 

The magnetic field penetration depth is a basic pa- 
rameter of superconductors, closely related to the super- 
fluid densitj* 2 -. Earlier on, the linear temperature depen- 
dence of the magnetic field penetration depth A(T) was 
observed for the cuprate superconductor YBa2Cu307_ y 
at low temperatures (T = 4K~20K)-, which first pro- 
vided a strong experimental support for the nodes in 
the d-wave SC gap function of cuprate superconduc- 



tors, then confirmed by the angle-resolved photoemission 
spectroscopy (ARPES) experiments^. Later, this lin- 
ear temperature dependence of the magnetic field pene- 
tration depth has been observed in different families of 
cuprate superconductors^—. However, at extremely low 
temperatures (T < 4K), the linear temperature depen- 
dence of the magnetic field penetration depth is modi- 
fied, and a nonlinearity emerges 1 ^—. Moreover, some 
indications of nonlocal effects giving rise to the non- 
linearity have been reported in the field dependence of 
the effective magnetic field penetration depth in cuprate 
superconductors! 4 . Furthermore, the doping dependence 
of the electromagnetic response in cuprate superconduc- 
tors has been studied in terms of the zero-temperature su- 
perfluid density. The superfluid density is proportional to 
the squared amplitude of the coherent macroscopic wave 
function describing the SC charge carriers, and there- 
fore it is an important physical quantity and can provide 
significant information about the SC state. In partic- 
ular, the superfluid density of cuprate superconductors 
in the underdoped regime vanishes more or less linearly 
with decreasing doping concentration^—. This in turn 
gives rise to the linear relation between the critical tem- 
perature T c and the superfluid density observed in the 
underdoped regime 1 ^. 

Theoretically, the electromagnetic response in cuprate 
superconductors has been extensively studied based 
on the the phenomenological Bardeen-Cooper-Schrieffer 
(BCS) formalism with the d-wave SC gap function-^—. 



2 



It has been shown 4 *^ that for a d-wave superconduc- 
tor in the local limit (£ <C A, where £ is the coher- 
ence length), the simple d-wave pairing state (assuming a 
tetragonal symmetry and ignoring the dispersion in the c- 
axis direction) gives the magnetic field penetration depth 
A(T) oc T/Ao, where Ao is the zero-temperature value 
of the d-wave gap amplitude. In particular, it has been 
argued that this linear temperature dependence of the 
magnetic field penetration depth is attributed to the ex- 
citation of quasiparticles out of the condensate at the 
nodes of the SC gap function. Furthermore, the fact£~— 
that this linear relation holds down to very low doping 
concentrations suggests that near the nodes these quasi- 
particle excitations are well described by a simple BCS- 
like formalism with the d-wave SC gap function, even for 
the doping concentration d — > 022. This is also consistent 
with the ARPES experiments^ 3 -. However, this depends 
sensitively on the quasiparticle scattering. In particu- 
lar, at extremely low temperatures, the coherence length 
may diverge at the nodes. This may imply that the local 
condition no longer holds, and the electromagnetic field 
varies significantly over the size of a Cooper pair. Conse- 
quently, the nonlocal effect emerges-- 3 and then plays an 
important role in the electromagnetic response of cuprate 
superconductors 1 ^—. It has been suggested 1 ^— that 
nonlocal effects can imply a crossover from the linear 
temperature dependence of the magnetic field penetra- 
tion depth at low temperatures to a nonlinear one in the 
extremely low temperature range. To the best of our 
knowledge, the electromagnetic response in cuprate su- 
perconductors has not been treated starting from a mi- 
croscopic SC theory, and no explicit calculations of the 
doping dependence of the superfluid density in the un- 
derdoped regime have been made so far. 

Recently, a kinetic energy driven SC mechanism has 
been developed 24 , where the charge carrier-spin interac- 
tion from the kinetic energy term induces a charge carrier 
pairing state with the d-wave symmetry by exchanging 
spin excitations. Then the electron Cooper pairs origi- 
nating from the charge carrier pairing state are due to 
charge-spin recombination, and their condensation re- 
veals the d-wave SC ground-state. In particular, this SC- 
state is controlled by both the SC gap function and the 
quasiparticle coherence, then the maximal SC transition 
temperature occurs around the optimal doping, and de- 
creases in both underdoped and overdoped regimes. The 
unique feature of this kinetic energy driven SC mecha- 
nism is that the pairing comes out from the kinetic en- 
ergy by exchanging spin excitations and is not driven by 
the magnetic superexchange interaction as in the reso- 
nant valence bond type theories 2 ^. Within the frame- 
work of the kinetic energy driven superconductivity, we 
have discussed the low energy electronic structur o 26 ' 27 of 
cuprate superconductors and the spin respons o 24 ' 28 , and 
qualitatively reproduced some main features of ARPES 
experiments^ and inelastic neutron scatterin g 29 ' 30 mea- 
surements on cuprate superconductors. 

The layered crystal structure gives rise to a strong 



anisotropy of cuprate superconductors, and it is possible 
to observe both in-plane and inter-plane electromagnetic 
responses. The former one is characterized by the ab- 
plane magnetic field penetration depth, whereas the lat- 
ter one is related to the magnetic field penetration in the 
c-axis direction. In this paper we concentrate on the in- 
plane electromagnetic response based on the kinetic en- 
ergy driven superconductivity and do not consider c-axis 
properties, which can be discussed, e.g., by taking into 
account hopping between adjacent copper-oxides layers 
within the tunneling Hamiltonian approach 22 . 

The paper is organized as follows. Within 
the framework of the kinetic energy driven d-wave 
superconductivity 2 ^, we discuss the electromagnetic re- 
sponse of cuprate superconductors in Section [TTl deriving 
the kernel of the linear response with a purely trans- 
verse vector potential. In Section Hm based on the spec- 
ular reflection mode l 31 ' 32 , we calculate the temperature 
and doping dependence of quantitative characteristics of 
the electromagnetic response, such as the local magnetic 
field profile, the magnetic field penetration depth, and 
the superfluid density. Our results show that the elec- 
tromagnetic response in cuprate superconductors can be 
understood within the framework of the kinetic energy 
driven d-wave SC mechanism in the presence of a weak 
external magnetic field. We conclude the paper with a 
brief summary in Section IIVI In Appendix [X] we present 
a method to generalize the approach in order to obtain 
gauge invariant results. 



II. ELECTROMAGNETIC RESPONSE IN 
CUPRATE SUPERCONDUCTORS 

A common feature of cuprate superconductors is the 
presence of two-dimensional Cu02 planes, and it is be- 
lieved that the unconventional physical properties of 
these systems are closely related to the doped Cu02 
plane ' . It has been argued that the essential physics of 
the doped Cu02 pian o 7 ' 25 is captured by the t-J model 
on a square lattice. However, for discussions of the elec- 
tromagnetic response in cuprate superconductors, the t- 
J model can be extended by including the exponential 
Peierls factors as, 

If/a la 

+ J^S,-S ;+f? , (1) 

if, 

where i) = ±x,±y, C\ a (C; CT ) is the electron creation 
(annihilation) operator, S; = (S[,Sf,S[) are spin op- 
erators, fi is the chemical potential, and the exponen- 
tial Peierls factors account for the coupling of electrons 
to the weak external magnetic field in terms of the vec- 
tor potential A(?) 33 ' 34 . This t-J model is subject to an 
important local constraint ^ a C\ a Ci a < 1 in order to 
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avoid the double occupancy. The strong electron cor- 
relation in the t-J model manifests itself by this local 
constraint^, which can be treated properly in analyti- 
cal calculations within the charge-spin separation (CSS) 
fermion-spin theory 3 ^, where the constrained electron op- 
erators are decoupled as Cif = h^Sf and Cq. = hj^Sf , 
with the spinful fermion operator hi„ = e~ l ® l ''hi repre- 
senting the charge degree of freedom together with some 
effects of spin configuration rearrangements due to the 
presence of the doped hole itself (charge carrier), while 
the spin operator Si represents the spin degree of free- 
dom. In particular, it has been shown that under the 
decoupling scheme, this CSS fermion-spin representation 
is a natural representation of the constrained electron 
defined in the Hilbert subspace without double electron 
occupancy^ 6 -. The advantage of this CSS fermion-spin 
approach is that the electron single occupancy local con- 
straint is satisfied in analytical calculations. Further- 
more, these charge carriers and spins are gauge invariant, 
and in this sense, they are real and can be interpreted as 
the physical excitations^ 6 .. In this CSS fermion-spin rep- 
resentation, the t-J model ([IJ can be expressed as, 

in 

+ h} +u hi±Sf S+ +fj ) 

- fl^hj^hicr + Jeff' X! S( ' S '+^' ( 2 ) 
la Irj 

where J c s = (1 — <5) 2 J, and S — {h\ a hi a ) = (hjhi) is the 
charge carrier doping concentration. As an important 
consequence, the kinetic energy term in the t-J model 
has been transferred as the charge carrier-spin interac- 
tion, which reflects that even the kinetic energy term in 
the t-J Hamiltonian has strong Coulomb contribution 
due to the restriction of no double occupancy of a given 
site. 

In the case of zero magnetic field, we^i have shown 
in terms of Eliashberg's strong coupling theory 3 -! that 
the charge carrier-spin interaction from the kinetic en- 
ergy term in the t-J model © induces a charge car- 
rier pairing state with the d-wave symmetry by exchang- 
ing spin excitations in the higher power of the charge 
carrier doping concentration 6, then the SC transition 
temperature is identical to the charge carrier pair tran- 
sition temperature. Moreover, it has been shown that 
this SC state is the conventional BCS-like with the d- 
wave symmetr} ^ 26 ! 27 , so that the basic BCS formalism 
with the d-wave SC gap function is still valid in quanti- 
tatively reproducing all main low energy features of the 
SC coherence of the quasiparticle peaks in cuprate su- 
perconductors, although the pairing mechanism is driven 
by the kinetic energy by exchanging spin excitations, and 
other exotic magnetic scatterin g 29 ' 30 is beyond the BCS 
formalism. Following the previous discussion s 24 ' 26 ' 27 , the 
full charge carrier diagonal and off-diagonal Green's func- 
tions in the SC state can be obtained explicitly in the 



Nambu representation as, 

rr°> /| ■ \ r, iu n T + £kT 3 - A hZ (k)ri 

<G(k, iw n ) = Z hF — -2 =2 , (3) 

where tq is the unit matrix, t\ and T3 are Pauli ma- 
trices, the renormalized charge carrier excitation spec- 
trum £k = ZhF^k, with the mean-field (MF) charge car- 
rier excitation spectrum £k = Ztx"fk — fJ-, the spin cor- 
relation function X = <S+S^>, 7k = (1/Z) £^ k ^, 
Z is the number of the nearest neighbor sites, the 
renormalized charge carrier d-wave pair gap function 
Ahz(k) = ZhFAh(k), where the effective charge car- 
rier d-wave pair gap function Ah(k) = Ah7 k <i ' > with 
7^ = (cosk x ~ cos/cj,)/2, and the charge carrier quasipar- 
ticle spectrum = \J^. + |Ahz(k)| 2 , while the charge 

carrier quasiparticle coherent weight Z^-p and the effec- 
tive charge carrier gap parameter Ah have been deter- 
mined self-consistently along with another seven quanti- 
ties and correlation function a 24 ' 26 ' 27 . Let us emphasize 
that the quasiparticle coherent weight renormalizing the 
physical quantities naturally emerges in our formalism 
(3), and then both the SC gap function and the quasi- 
particle coherence control the SC state. Therefore in our 
approach there is no need to introduce any phenomeno- 
logical charge renormalization factors in order to describe 
the electromagnetic response^ 2 .. 

Now we turn to the discussion of the electromagnetic 
response in the kinetic energy driven cuprate supercon- 
ductors. The weak external magnetic field applied to 
the system usually represents a weak perturbation, but 
the induced field generated by supercurrents cancels this 
weak external field over most of the volume of the sam- 
ple. Consequently, the net field acts only very near the 
surface on a scale of the magnetic field penetration depth 
and so it can be treated as a weak perturbation on the 
system as a whole. Therefore the electromagnetic re- 
sponse can be successfully studied within the linear re- 
sponse approac h 38 ' 39 , where the averaged value J of the 
induced microscopic screening current j in the presence 
of the vector potential A is found as, 

3 

Jp, (q, w) = - K ^ (q> u ) A v (q> w ) ) ( 4 ) 

!/=l 

where is the kernel of the response function and 
the Greek indices label the axes of the Cartesian coordi- 
nate system. Recall that, as always in the linear response 
method, the thermal average of the supercurrent is calcu- 
lated with the unperturbed Hamiltonian, i.e. for A = 
in Eq. ((2]). Let us also note that the relation ((4]), which is 
local in the reciprocal space, in general implies a nonlocal 
response in the coordinate space. 

The kernel, which plays a central role in the description 
of the electromagnetic response, and once known allows 
us to calculate quantitative characteristics of the electro- 
magnetic response, can be separated into two parts: 

K^fau) = K${q,u) + K$(a,u), (5) 
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a diamagnetic part Kpt) and a paramagnetic one K^J . 
The evaluation of the diamagnetic contribution usually 
poses no difficulties since it is known almost immediately 
from the form of the diamagnetic current operator: it 
turns out to be diagonal and proportional to the average 
kinetic term. However, the paramagnetic part can only 
be calculated approximately since it involves evaluation 
of a retarded current-current correlation function (polar- 
ization bubble). As the retarded function is inconvenient 
for perturbation analysis one usually proceeds with the 
corresponding imaginary-time-ordered Matsubara func- 
tion, 



P„ 



-<r T {i(P)(q,T)jW(-q,0)}), 



(6) 



where the paramagnetic current operator is defined in the 
imaginary time t Heisenberg picture. Hence, the main 
problem is reduced to the evaluation of a retarded current 
commutator for the unperturbed system. The retarded 
current-current correlation function is then obtained in a 
standard way from the imaginary time Fourier transform 
P Mt/ (q, iu> n ) of the Matsubara function ([6]) by analytic 
continuation to real frequencies 3 ^. 

The vector potential A (then the weak external mag- 
netic field B = rotA) has been coupled to the elec- 
trons, which are now represented by Cif = hJ^Sf and 
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(j. = hnSi~ in the CSS fermion-spin representation. 
However, in the CSS framework, the vector potential A 
is coupled to h, , while the corresponding weak external 
magnetic field B = rotA is coupled to S; by including 
the Zeeman terrr42 in the Hamiltonian (TTJ). For cuprate 
superconductors, the upper critical magnetic field is 50 
Tesla or greater around the optimal doping. In this pa- 
per, we mainly focus on the case where the applied ex- 
ternal magnetic field B < 10 mT is much less than the 
upper critical magnetic field. In this case, the Zeeman 
termiS in the Hamiltonian fl} has been dropped, and 

then the electron current operator — + can 
be obtained by differentiating the Hamiltonian ([2]) with 
respect to the vector potential as, 



.•(d) 



,-(p) 



^EK^^ + ^^m)^^), (7a) 
i\et < 

la 



la 



(7b) 



being the diamagnetic and paramagnetic contributions, 
respectively. 

Since the diamagnetic current is explicitly proportional 
to the vector potential, it is straightforward to find the 



diamagnetic part of the response kernel as, 



4 d , } (q) 



Z hFX e 2 t 1 



h 2 



N 

k 



<L„ cos k u 



, 6c , , 

1 — — — tanh — - — 

£-hk 2 



(8) 



h 2 ^ 

The paramagnetic part of the response kernel is more 
complicated to calculate, as it involves evaluation of the 
current-current correlation function (J6)). In particular, if 
we want to keep the theory gauge invariant, it is crucial 
to approximate the correlation function in a way main- 
taining local charge conservatioit^ i 34 ' 39 i 41 . Since in the 
following calculations we will work with a fixed gauge 
of the vector potential, we postpone the detailed discus- 
sion of this problem until Appendix |XJ Starting with 
the paramagnetic current operator (|7bp . we can rewrite 



its Fourier transform in the notation of Nambu ^fields 
" - (u ' [ ' ■ ^ and* k+q = f/ik+ q t^L k _ q J 



l k+qf> 



as 



4 p) (q) 



N Z—i k 



~ e l 2 sm 
n 



2 



To 



(9) 

For the purpose of the discussion addressing the gauge 
invariance problem, presented in Appendix 13 it is con- 
venient to find the charge density in the Nambu no- 
tation as well. Within the CSS fermion-spin scheme, 
we first find p(q) ps -(e/2JV)^ k (<5 q ,o - ^t^k+qt - 
/i^/jk+q^)- Then the paramagnetic four-current oper- 
ator can be represented in the Nambu form as (q) = 
J2 ^k7ft(k + q, k)<3> k+q , where the bare current vertex, 



7 M (k + q,k) 



Xet 



e l2 £ sin (kfj, + ^) r for ^ ^ 



■ T3 for n = 0. 

(10) 

It is necessary to be aware that we are calculating the 
polarization bubble with the paramagnetic current op- 
erator ([9]), i.e., bare current vertices (fl0|) . but charge 
carrier Green functions. Consequently, as in this sce- 
nario we do not take into account longitudinal excita- 
tions properly^ 4 -, the obtained results are valid only in 
the gauge, where the vector potential is purely trans- 
verse, e.g. in the Coulomb gauge. In this case, we can 
obtain the correlation function © in the Matsubara rep- 
resentation as, 



J 



P (tl!/ (q,iw„) = — 



Xet 
h 



1 



K9m-9"-)_ > sm 
N ^ 

k 



+ D sin + I £Tr[G(k 



+ q, iu n + iu m )G(k, iv m )] . (11) 
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Restricting the discussion to the static limit (ui ~ 0) 
and completing the summation over Matsubara frequen- 
cies, we obtain the bare vertex current-current correla- 



tion function, and hence the paramagnetic part of the 
response kernel as, 



42 M) 



\ h 2 J 
l 

Ehk + Shk+q 
1 

£-hk — -Elik+q 



N ^ 



sin ( k 



sin 



(*» + f) 



_ Ck+qCk + A hZ (k + q)Ahz(k) 

-Ehk-Elik+q 

1 | £k+q£k + A hZ (k + q)A hZ (k) 



[1 - n F (i?hk) - ^(Shk+q)] 
[nF(Shk+ q ) - n F (E h k)} 



(12) 



Note that in the long wavelength limit, when |q| — > 0, 
the former term in Eq. (|12p vanishes, and the latter 

turns into -2 (xetZ hF /Nti 2 ) £ sinfc M sin k v np(£ ; hk)[l — 

k 

np(-Ehk)], which is equal to zero in the zero-temperature 
limit. Hence, in this case, the long wavelength electro- 
magnetic response at zero temperature is determined by 
the diamagnetic part of the kernel only. 



III. QUANTITATIVE CHARACTERISTICS 

The way the system reacts to a weak electromag- 
netic stimulus is entirely described by the linear response 
kernel, which is calculated within a microscopic model. 
Once the kernel is known, the effect of a weak external 
magnetic field can be quantitatively characterized by ex- 
perimentally measurable quantities such as the magnetic 
field penetration depth and the local magnetic field pro- 
file. Technically, we need to combine one of the Maxwell 
equations with the relation ((4]) describing the response 
of the system and solve them together for the vector 
potential. This is the step in which a particular gauge 
of the vector potential — usually implied by the geome- 
try of the system — is set. However, the kernel function 
derived within the linear response theory describes the 
response of an infinite system. In order to take into ac- 
count the confined geometry of cuprate superconductors 
it is necessary to introduce a surface being the boundary 
between the environment and the sample. This can be 
done within the standard specular reflection mode l 31 i 32 
with a two-dimensional geometry of the SC plane, in the 
configuration with external magnetic field perpendicular 
to the ab plane, as shown in Fig. [TJ In the present paper 
we study magnetic field penetration effects within the ab 
plane only, so our primary goal is to find and discuss the 
magnetic field in-plane penetration depth. 

In order to simulate an external magnetic field at 
the surface of a two-dimensional sample, we introduce 
an external current sheet J y , ox t(x) = — 2BoS(x)/ /iq at 
the edge x = 0, where /xq is the magnetic permeabil- 



external field 
B = (0, 0, B ) 




total current 

J M = (0, J vM +J„ 



- 0),-'' 



X - penetration depth 

FIG. 1: Geometry of the specular reflection model. The cur- 
rent Jcxt simulates external magnetic field at the edge of the 
sample (x — 0), whereas the induced supercurrent Ji nt is the 
(linear) reaction of the system. 



ity and Bq is the amplitude of the weak external mag- 
netic field at the surface (x = 0). From the Maxwell 
equation for the curl of the local magnetic field rot h = 
Mo (Jint + Jext) = Mo Jint + [0, -2B 5(x), 0] and the fact, 
that the induced supercurrent J; n t flows along the y axis, 
we can state that the local magnetic field is of the form 
h(r) = [0, 0, h z (x)}. In order to discuss the magnetic field 
penetration effect, spatial dependence of the local mag- 
netic field has to be found. Let us begin with the iden- 
tity rot rot A = graddivA — V 2 A and choose the vec- 
tor potential as A(r) = [0, A y (x), 0] setting the Coulomb 
gauge. In this case, q 2 x A y {q L ) = fi [J y , in t(q) + J y ,cxt(q)] , 
because the vector potential has only non-zero y compo- 
nent. Finally, including the form of the external current, 
the linear relation Q between the induced supercurrent 
and the vector potential J y ,i n t(q) = —K yy (q)A y (q), and 
solving for the vector potential we obtain, 



A y (q) = -8ir 2 B 



K%) 5 {lz) 



(13) 



Since the vector potential has only the y component, 
the only non-zero component of the local magnetic field 
h = rot A is that along the z axis and h z (q) = iq x A y (q). 
Substituting the derived form of the vector potential (|13l) , 
and taking the inverse Fourier transform, the local mag- 
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netic field profile can be obtained as 
dq x 



n z [x) = — 

7T 



q x smq x x 



^oK yy {q Xl O,0) + 



(14) 



Local magnetic field profiles can be measured experimen- 
tally, e.g. using the muon-spin rotation techniqu e) 12 ^ 13 , 
providing an important tool to investigate the details of 
magnetic field screening inside the sample. In cuprate 
superconductors the screening is found to be of exponen- 
tial characte r 12 ! 13 , in support of a local (London-type) 
nature of the electrodynamics 2 -. For the convenience of 
the following discussions, we introduce a characteristic 
length scale ao = \J h 2 a / ^e 2 J . Using a reasonably 
estimative value of J/fce ~ 1000K and a « 0.383nm, 
which is the lattice parameter for the cuprate super- 
conductor YBa2Cu307_j,, we obtain ao ~ 97.8nm. In 
Fig. [2] we plot the local magnetic field profile (fT4)) as a 
function of the distance from the surface at temperature 
T = 0.02J for the doping concentration S — 0.150 (solid 
line), S = 0.147 (dashed line), and S = 0.144 (dotted 
line) with parameter t = 2.5 J. For comparison, the cor- 
responding experimental result of the local magnetic 
field profiles for the high quality YBa2Cu307_ a sample 
is also shown in Fig. [2] (inset, bottom). If a weak ex- 
ternal field Bo » 10 mT is applied to the system just 
as it has been done in the experimental measurement^, 
then the experimental result 13 for YBa2Cu307_ a is well 
reproduced. In particular, our theoretical results per- 
fectly follow an exponential law as expected for the local 
electrodynamic response. 



m 




160 



FIG. 2: The local magnetic field profile as a function of the 
distance from the surface at temperature T — 0.02J for dop- 
ing concentration S = 0.150 (solid line), 8 = 0.147 (dashed 
line), and 8 = 0.144 (dotted line) with parameter t — 2.5 J. 
Inset (top): zoom into the intermediate range of the local 
magnetic field profile. Inset (bottom): the corresponding ex- 
perimental result for YBa2Cu307_ a taken from Ref. fl3l . 

The above obtained local magnetic field profile h z (x) 
allows us to determine the magnetic field in-plane pene- 
tration depth A(T) in a straightforward way. According 
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FIG. 3: Temperature dependence of the magnetic field in- 
plane penetration depth AA(T) for the doping concentration 
5 = 0.150 (solid line), 8 = 0.149 (dashed line), and S = 0.148 
(dotted line) with parameter t/J — 2.5. Inset: the corre- 
sponding experimental data for YBa2Cu307_j / taken from 
Ref. i. 



to the definition \(T) = B 1 J °° h z {x) da;, the magnetic 
field in-plane penetration depth can be evaluated as, 

oo 

dgz 



KT) = ~ [ 

7T J 



fi o K y y(q x ,0, 0) 



(15) 



In this case, we obtain the zero-temperature magnetic 
field in-plane penetration depth A(0) « 380. 8nm for 
the doping concentration <5 = 0.150 with parameter 
t/J = 2.5. This anticipated value is very close to 
the values of the magnetic field in-plane penetration 
depth A ~ 156nm ~ 400nm observed in different fam- 
ilies of cuprate superconductor a 12 ' 17 i 42 . Furthermore, 
AA(T) = A(T) — A(0) as a function of temperature T for 
the doping concentration S = 0.150 (solid line), S = 0.149 
(dashed line), and 5 = 0.148 (dotted line) with parameter 
t/J = 2.5 is plotted in Fig. [3]in comparison with the cor- 
responding experimental results^ of YBa2Cu307_ y (in- 
set). Our theoretical results show linear characteristics 
of the magnetic field in-plane penetration depth AA(T), 
except for extremely low temperatures where a strong 
deviation from the linear characteristics (a nonlinear ef- 
fect) appears. In particular, this crossover from the linear 
temperature dependence in the low temperature regime 
into the nonlinear one at extremely low temperatures 
is observed experimentally in nominally clean crystals 
of cuprate superconductors^ - — . Apparently, there is 
a substantial difference between theory and experiment, 
namely, the value of the difference between A(T) and 
A(0) calculated theoretically is much smaller than the 
corresponding value measured in the experiment. How- 
ever, upon a closer examination one can see immediately 
that the main difference is due to fact that the calcu- 
lated A(T) increases slowly with temperature. As for a 
qualitative discussion in this paper, the overall tendency 
seen in the theoretical result is consistent with that in the 
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experiment^. In cuprate superconductors, the values of J 
and t are believed to vary somewhat from compound to 
compound 7 -. Therefore the quantitative agreement can 
be reached by adjustments of theory's parameters t and 
J, or by introducing the next neighbor hopping t' . 
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FIG. 4: Doping dependence of the zero-temperature in- 
plane superfluid density in the underdoped regime with 
t/J — 2.5. Inset: the corresponding experimental result for 
YBa2Cu307- H taken from Ref. 0. 

A quantity which is closely related to the magnetic 
field in-plane penetration depth A(T) is the in-plane su- 
perfluid density p s (T) = A _2 (T). For a better under- 
standing of the physical properties of cuprate supercon- 
ductors, we have calculated the doping dependence of 
the zero-temperature in-plane superfluid density p s (0) 
in the underdoped regime. The result for parameter 
tj J = 2.5 is plotted in Fig. [4] in comparison with the cor- 
responding experimental data 16 for YBa2Cu307_ a (in- 
set). It is shown that the in-plane superfluid density 
p s (0) in the underdoped regime vanishes more or less lin- 
early with decreasing doping concentration <5, in qual- 
itative agreement with experimental results of cuprate 
superconductors^—. This result also is a natural conse- 
quence of the linear doping dependence of the SC transi- 
tion temperature T c oc 6 in the underdoped regime in the 
framework of the kinetic energy driven SC mechanism 24 , 
where the SC transition temperature T c is set by the 
charge carrier doping concentration, and then the den- 
sity of the charge carriers directly determines the in-plane 
superfluid density in the underdoped regime. 

The appearance of the nonlinearity in the tempera- 
ture dependence of the magnetic field in-plane penetra- 
tion depth in cuprate superconductors at extremely low 
temperatures, as shown in Fig. [3l can be attributed to 
the nonlocal effects, which in the case of a pure d-wave 
cuprate superconductor with nodes in the gap become 
significant for the electromagnetic response^—. In gen- 
eral, the relation between the supercurrent and the vec- 
tor potential ([J| is nonlocal in the coordinate space due 
to the finite size of charge carrier Cooper pairs. In the 
framework of the kinetic energy driven d-wave SC mech- 
anism, the size of charge carrier pairs in the clean limit is 
of the order of the coherence length £(k) = fa; F /7rAh(k), 



where v-p — ft~ 1 9^k/9k|/ CF is the charge carrier velocity 
at the Fermi surface, and therefore the size of charge car- 
rier pairs is momentum dependent. Although the weak 
external magnetic field decays exponentially on the scale 
of the magnetic field in-plane penetration length A(T), 
any nonlocal contributions to measurable quantities are 
of the order of k~ 2 , where k, known as the Ginzburg- 
Landau parameter, is the ratio of the magnetic field 
in-plane penetration depth A and the coherence length 
£. However, in the d-wave cuprate superconductors, 
the characteristic feature is the existence of four nodal 
points [±7r/2, ±7r/2] in the Brillouin zone, where the 
charge carrier gap function vanishes Ah(k)|^ 7r / 2 ,±ir/2] — 
Ah(cosfca; — cos/cj,)/2|[± 7r / 2 .±7r/2] =0. As a consequence, 
the coherence length £(k) diverges around the nodes. 
In particular, at extremely low temperatures, the quasi- 
particles selectively populate the nodal region, and the 
major contribution to measurable quantities comes from 
these quasiparticles. In this case, the Ginzburg-Landau 
ratio Ai(k) around the nodal region is no longer large 
enough for the system to belong to the class of type- 
II superconductors, and the condition of the local limit 
is not fulfilled^. On contrary, the system falls then 
into the extreme nonlocal limit, and therefore the non- 
linear characteristic in the temperature dependence of 
the magnetic field in-plane penetration depth can be ob- 
served experimentally in cuprate superconductors at suf- 
ficiently low temperature o 3 ' 12 " — . However, with increas- 
ing temperature, the quasiparticles around the nodal re- 
gion become excited out of the condensate, and the non- 
local effect fades away. In this case, the momentum de- 
pendent coherence length can be replaced approxi- 
mately with the isotropic one £o = Hy-p/irA^. Then the 
Ginzburg-Landau parameter kq ps A(0)/(o ~ 180, and 
the condition for the local limit is satisfied. This antici- 
pated value of the Ginzburg-Landau parameter kq ss 180 
is not too far from the range kq w 150 ~ 400 esti- 
mated experimentally for different families of cuprate 
superconductor a 12 ! 17 ' 42 . Consequently, the cuprate su- 
perconductors at moderately low temperatures turn out 
to be type-II superconductors, where nonlocal effects are 
negligible, the electrodynamics is purely local and the 
magnetic field decays exponentially over a length of the 
order of a few hundreds nm. In this local limit, the 
pure d-wave pairing state in the kinetic energy driven 
SC mechanism gives the magnetic field penetration depth 
A(T) cx T— This is why the linear temperature depen- 
dence of the magnetic field in-plane penetration depth 
A(T) is observed experimentally 3 ^^— in cuprate super- 
conductors at moderately low temperatures. 

Finally, we have to note that a deviation from the lin- 
ear Uemura relation between the in-plane superfluid den- 
sity p s (0) and doping concentration 5 has been observed 
recently in the underdoped regim o 11 ' 1 ^ 43 . This deviation 
from the linear Uemura relation suggests a sublinear de- 
pendence of the critical temperature T c and the super- 
fluid density p s (0), since T c must fall to zero when p s (0) 
doe o 11 ' 43 . The parent compound of doped cuprate super- 
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conductors is a Mott insulator with an antiferromagnetic 
long-range order and superconductivity occurs when the 
antiferromagnetic long-range order state is suppressed by 
doped charge carriers. Since these doped charge carriers 
in cuprate superconductors are induced by the replace- 
ment of some ions by other ones with different valences, 
or the addition of excess oxygens in the block layer, there- 
fore, in principle, all cuprate superconductors have nat- 
ural impurities 7 . Therefore the impurities play an im- 
portant role in the electromagnetic response and lead to 
some subtle differences in the electromagnetic response 
for different families of cuprate superconductors — . In this 
case, the impurity effect on the SC state of cuprate super- 
conductors is also a possible source for the deviation from 
the linear Uemura relation. In this context we^ have 
discussed the effect of the extended impurity scatterers 
on the quasiparticle transport of cuprate superconduc- 
tors in the SC state based on the nodal approximation 
of the quasiparticle excitations and scattering processes, 
and predicted that in contrast with the dome shape of 
the doping dependent SC gap parameter, the minimum 
of the microwave conductivity occurs around the opti- 
mal doping, and then increases in both underdoped and 
overdoped regimes. However, in this paper we are pri- 
marily interested in exploring the general notion of the 
electromagnetic response in cuprate superconductors in 
the SC state. The qualitative agreement between the 
present theoretical results in the clean limit and experi- 
mental data for different families of cuprate superconduc- 
tors provides an important confirmation of the nature of 
the SC phase of cuprate superconductors as a d-wave 
BCS-like SC state within the kinetic energy driven SC 
mechanism. 



IV. CONCLUSIONS 

In this paper we have discussed the electromagnetic 
response in cuprate superconductors within the frame- 
work of kinetic energy driven d-wave superconductivity. 
Following the linear response theory and taking into ac- 
count the two-dimensional geometry of cuprate super- 
conductors within the specular reflection model, we have 
reproduced some main features of the electromagnetic 
response experiments on cuprate superconductors, in- 
cluding the exponential local magnetic field profile, the 
linear temperature dependence of the in-plane penetra- 
tion depth in the low temperature range and its nonlin- 
ear temperature dependence at extremely low temper- 
atures. Moreover, the linear doping dependence of the 
zero-temperature in-plane superfluid density in the un- 
derdoped regime has been reproduced. In particular, we 
have clearly identified the limitations of the used approx- 
imations, especially with respect to the problem of gauge 
invariance. Furthermore, we have proposed a method to 
generalize the discussions in order to make them inde- 
pendent of a particular choice of the vector potential. 
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Appendix A: Gauge-invariant electromagnetic 
response 

It is well known that gauge invariance is a direct con- 
sequence of local charge conservation^ 3 ^, which is math- 
ematically expressed by the charge density-current conti- 
nuity equation or its Green function analogue called the 
generalized Ward identity (GWI^^^. 

3 

-2iV^ (?M r Al (fc+g,fc) -r 3 G- 1 (fc)-G- 1 (fc+g)r3. (Al) 

Here T M is a dressed version of the four-current vertex 
function, and the four- vector notation q = (q, go = 
along with the metric (1, 1, 1, —1) has been introduced. 

Since the local charge conservation requirement is quite 
universal and fundamental, it should be inherent to any 
theory of the electromagnetic response which is expected 
to be gauge invariant. The purpose of this appendix is 
to propose — within the framework of the kinetic energy 
driven superconductivity — a method to dress the current 
vertex in a way, which does not violate the GWI. Once 
such a method is found, the bare polarization bubble 
([6]) can be replaced with its dressed version presented in 
Fig. [5J and the resulting kernel of the response function 
will provide correct results for any gauge of the vector 
potential. 



G 




FIG. 5: Dressed polarization bubble (Nambu notation). Here 
both the Green function and the current vertex are dressed 
with the pairing interaction due to the spin bubble. 

In the first step we will note that 

3 

-2N J2 <7m7m (k + q,k) = ^G^-^-G^-^k + q)^ 

fj,=0 

(A2) 
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i.e. that the GWI for the bare current vertex is 
satisfied with the MF charge carrier Green function 
GW(A) = [(iuj n ) 2 - Ck] _1 (^«r + 6cT 3 ). Substitut- 
ing the MF charge carrier Green function, the rhs of 
Eq. (HH) turns into TgG^^^fc) - G^" 1 ^ + q)r 3 = 
(Ck+q — £k) To — Moreover, in the long wave- 

length limit, after including the explicit form of the MF 
charge carrier dispersion relation found within the frame- 
work of the kinetic energy driven superconductivity^, 
it further simplifies to r 3 G(°) -1 (A:) - G (0)_1 (fc + q)r 3 « 
[—2tx(q x sinfca; + q y smk y )} tq — qoT 3 . Now, recalling the 



form of the bare vertex (fTOf . we can notice that 
in the long wavelength limit the scalar product on 
the left-hand side of Eq. (|A2|) — <Zo7o + q7 = 
(27V) -1 (r 3 <7o — 7oVk£k • q), which proves the equality 

It is well known that in order to obtain a dressed vertex 
function, which does not violate the GWI, a ladder-type 
approximation can be adapted^ ' 34 i 39 . The nature of the 
pairing mechanism- 4,26 , which originates from the spin 
bubble, suggests a ladder-like approximation of the form, 



1 1 



T^{k + q,k) = ltl {k + q,k) + — -J2r3G(k+ P + q)T^{k+ P + q,k+p)G(k+p)T 3 



N H A p+ P '+k n (p> p'; iPm), 



(A3) 



which is graphically presented in Fig. [6] 




FIG. 6: Ladder-type approximation for the dressed vertex. 



In order to prove that the approximation (| A3|) for 
the dressed vertex in fact implies a gauge invariant de- 
scription of the electromagnetic response, it is neces- 
sary and sufficient to check whether it does not vio- 
late the GWI (|A1|) . In order to prove it, we insert the 
dressed vertex function (|A3I) into the left-hand side of Eq. 
(|Alj) and use the identity — 2AT^ q^T^(s + q,s) = 
T3G _1 (s) — G _1 (s + q)r 3 to obtain 



3 3 1 1 / 1 \ 

^q^Y^k + q.k) = ^2<lnln(k + q,k) + — [ ) [T3<G(k + p + q)r 3 -T 3 



M=0 



li=0 



2N 



■p)vi\ 



A P+p'+k n (P ; P'; *Pm). 

p' 



(A4) 



In the long wavelength limit we use the approxima- 
tion A p+p , +k w A p+p/+k+q . Then we can simplify Eq. 

<EH> into £ 3 =0 q^k + q 7 k)^ £ 3 =Q q^{k + q,k)- 
[£(& + q)r 3 - T 3 E(fc)] /2N. Using the fact that the free 
vertex satisfies the GWI with the MF Green function, as 
stated in Eq. (|A2[) . and arranging the terms with respect 
to the Pauli matrices, we have 

3 

-2NY,QM k + 9,k) « T 3 [G^-\k)-J:(k)} 

(1=0 

- [G^- 1 (k + q)-E(k + q)}r 3 . 



Hence, identifying the terms in the square brackets as 
dressed charge carrier Green functions, we eventually ob- 
tain the GWI (|A"T1) . what proves that the ladder-type ap- 
proximation (|A3[) for the vertex function in the dressed 
polarization bubble in Fig. [5]is consistent with the GWI. 
Consequently, the kernel of the linear response calculated 
with the dressed polarization bubble is gauge invariant. 
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